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PREFACE 


This report is one of three reports describing results of investigations 
for the "Design, Fabrication, Exploratory Operation and Operational Manual of 
a Hydraulic Analog Computer . " The investigation has been sponsored by the 
Arctic Construction and Frost Effects Laboratory, New England Division, 

Corps of Engineers, U. S„ Army, under contract No. DA- 19 -016- Eng- 2743 with 
the Massachusetts Institute of Technology through its Division of Industrial Co- 
operation under D. I. C„ Project No. 5-7155. 

The entire investigation is directed toward improvement of techniques 
for predicting subsurface temperatures, particularly with reference to frost 
and thaw penetration below airfield pavements. This report presents the re- 
sults of studies related to frost penetration in multilayer soil profiles. Two 
other reports are (a) Design and Operation of an Hydraulic Analog Computer for 
Studies of Freezing and Thawing of Soils, and (b) Heat Transfer at the Air- 
Ground Interface with Special Reference to Airfield Pavements. 
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SYNOPSIS 


This report presents the results of an investigation relative to the pre- 
diction of the depth of frost penetration in multilayer soil profiles. An adapta- 
tion of the Modified Berggren formula is presented and computed results are 
compared with " exact” solutions for the depth of freezing as determined by the 
hydraulic analog computer developed under the contract. 

Nine cases were investigated, representing typical soil profiles below 
bituminous concrete airfield pavements . In seven cases, the pavement surface 
temperature was assumed to drop suddenly to a given value below freezing 
where it remained constant. In the remaining two cases, the surface tempera- 
ture was assumed to vary sinusoidally over an annual period. 

Results of these tests indicate that the proposed adaptation of the Modi- 
fied Berggren formula underpredicts the frost depth by approximately 3 percent. 
However, this error is believed to be insignificant compared to those errors in 
the predicted depth of frost penetration in any practical case arising from un- 
certainties in the thermal properties of the soil, particularly the thermal con- 
ductivity, and from unknown pavement surface temperatures. These factors 
present the most fruitful areas for future research on the depth of frost pene- 
tration prediction . 


xi 



I. INTRODUCTION 


1-01. CONTRACT 1 

On 2 October, 1953, the Massachusetts Institute of Technology through 
its Division of Industrial Cooperation entered on a cost-type research contract 
with the Arctic Construction and Frost Effects Laboratory, New England Divi- 
sion, Corps of Engineers, U. S. Army. This contract, No. DA-19-016-Eng- 
2743, is entitled "Design, Fabrication, Exploratory Operation and Operational 
Manual of a Hydraulic Analog Computer". This is one of three reports cover- 
ing investigations at M.I.T. during the first two years of this contract.' 

Specifically, this report describes results of investigations included 
within the following paragraphs of Appendix "A" of the contract: 

"3. After completing the work enumerated in Paragraph 1 and 2 above, 
the contractor shall perform comprehensive studies of the oper- 
ational use of the apparatus for several or all of the following in- 
vestigations: (a) Make a detailed study of the effect of variables, 
such as stratification, frost heave, and surface temperature fluc- 
tuations on the depth of freezing and thawing; compare the results 
of these studies with various formulas and averaging techniques 
which have been composed to handle these complicated situations." 

1-02. SCOPE 

The investigation described herein is directed primarily at evaluating 
the effectiveness of an adaptation of the Modified Berggren formula for predict- 
ing the depth of frost penetration in multilayer soil profiles. It is assumed, a 
priori, that the thermal properties of the soil and the boundary temperature 
conditions are known. The "exact" hydraulic analog solution for a given soil 
profile is compared with the predicted solution for these assumptions. 

Comparisons are made for both a step change in surface temperature 
and for a surface temperature varying sinusoidally to give a total of nine solu- 
tions for six different soil profiles, all having a 3 -inch bituminous concrete 
surface . 
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1-03, M.I.T, PERSONNEL 

Investigations under this contract have been carried out under the super- 
vision of Harl P. Aldrich, Jr., Assistant Professor of Soil Mechanics, Depart- 
ment of Civil and Sanitary Engineering. 

Mr. Rasmus S. Nordal, formally half-time Research Assistant in Civil 
Engineering, carried out the majority of the investigation described herein. 

Dr. Aldrich and Mr, Nordal are responsible for the technical writeup. 



II. GENERAL 


Freezing and thawing below highway and airfield pavements presents 
problems of utmost importance to their design and maintenance in Northern 
climates. These problems have assumed increased importance as wheel load- 
ings, traffic frequency, and cost of pavement construction and maintenance 
have increased. Furthermore, thermal problems in soil are receiving addi- 
tional attention as a result of construction activity in permafrost areas. 

It is generally recognized that three requirements must be met simul- 
taneously for significant ice -segregation and frost heave to occur in a foundation 
soil: (1) below freezing temperatures must penetrate into the ground, (2) the 
soil must meet certain grain- size requirements to be "frost- susceptible" and 
(3) a source of water must be available. This report deals with the first re- 
quirement only, namely the prediction of the rate and maximum depth of frost 
penetration through stratified (multilayered) soil profiles beneath bituminous 
concrete pavements. 

The problem of predicting the depth of frost penetration has been studied 
by numerous investigators during the last 30 years especially. All practical 
analytic solutions are derived for uniform soil. The most simple of these is the 
well-known Stefan formula: 



which may be written in terms of the freezing index F, expressed in degree 
days, as follows: 



( 2 - 2 ) 


This expression neglects the volumetric heat of the soil and therefore generally 
predicts depths of frost penetration which are too large. 

The Modified Berggren formula (1) (2) has been developed to account 


^Numbers in parenthesis refer to references in the Bibliography. 
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for volumetric heat as well as latent heat and may be written: 

./^v 

x = x V/^JL (2r3) 

where: 

X; (lambda); dimensionless correction coefficient which is given in 
Figure 14 as a function of two important dimensionless parameters, 
ot (alpha) and jjl (mu): 

v 

Thermal ratio, ca = — 

v 

s 

Fusion parameter, jo. = j- v g 

Introducing F as before to yield an equation for practical use, we obtain: 

X = X (2-4) 

Numerous other expressions (4) (5) (8) have been presented for frost penetration 
in uniform soil under one- dimensional conditions. 

No exact analytic solutions to the depth of frost penetration in multilayer 
soil profiles exist. Therefore, for this nonuniform case which occurs typically 
below highway and airfield pavements, semi -empirical computation techniques 
must be used. Since no exact solution for the multilayer case exists, there is 
no possibility of verifying, in a simple way, the reliability of the approximate 
procedure itself. 

It is true that a numerical approach based on a finite difference approxi- 
mation to the fundamental differential equations involved, will lead after con- 
siderable effort to an "exact' r solution. However, with the construction of the 
hydraulic analog computer (6) a relatively simple tool becomes available for 
checking semi- einpirical computational procedures. 

Thus, we seek to compare a series of solutions obtained on the hydraulic 
analog computer with results of approximate mathematical procedures for a 
variety of multilayer soil profiles which may occur below pavements. The first 
part of the report, however, presents a study of the accuracy of the analog com- 
puter solution itself. Results are compared with known analytic solutions for a 
step change in surface temperature and with results of numerical analyses car- 
ried out on the IBM -701 computer for a surface temperature varying sinusoidally. 



III. HYDRAULIC ANALOG COMPUTER 

3-01 o GENERAL 

/ 

The hydraulic analog computer developed under this contract is de- 
scribed in detail in Reference (6), Only a brief description of the computer will 
be presented here after which its operation relative to this investigation is given. 

The computer consists of a series of glass tube standpipes with expan- 
sion reservoirs about mid- height. Each finite layer of soil in a given subsurface 
profile is represented by one of these standpipes . Water level in a standpipe 
corresponds to temperature for the soil layer it represents. Flow of heat with- 
in the soil is represented in the computing device by the flow of water. The 
water level versus stored water characteristic of each standpipe and expansion 
reservoir is analogous, then, to the thermal energy versus temperature charac- 
teristic of the soil . 

To represent resistance to heat flow between adjacent soil layers, a re- 
sistance to fluid flow between standpipes must be provided. In the hydraulic 
analog computer, laminar flow in a hydraulic resistance is maintained in open- 
ings of capillary size, which can be continuously varied in length to yield the 
desired resistance . 

Before an actual problem can be solved on the computer, the soil profile 
and its thermal properties must be known or assumed. The computer can then 
be set up with the appropriate standpipes, latent heat reservoir openings, and 
hydraulic resistances. Ground surface temperature, represented by a water 
reservoir, can be programmed automatically to follow actual or assumed varia- 
tions in surface temperature with time by means of a moving paper template and 
servo-motor system. 

The programming and scaling procedure as well as photographs and 
drawings for the hydraulic analog computer are given in Reference (6), 

3 -02 . CHECK SOLUTIONS, UNIFORM SOIL 

Because the hydraulic analog computer solves a thermal problem by 
’Tumping” the continuous subsurface soil profile into a finite number of elements, 
the solution is not exact in the strictest sense. This "lumping” error depends, 
among other things, on the size of each lump and on the magnitude of thermal 
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gradients. These lumping errors were investigated by Scott (9). Since other 
minor errors are inevitable in the analog system, it becomes important to check 
its performance with known analytic solutions in order to gain confidence in the 
computer results. 

A series of four analog solutions were run which demonstrate that com- 
puter errors are negligible as compared to the precision of the physical data 

available for a practical case. A summary of these tests, numbered 1 through 

* 

4, is given in Table I and Figures 1 through 4. All tests were run assuming a 
homogeneous soil profile. Tests 1 and 2 were run under a step change in surface 
temperature while tests 3 and 4 are corresponding tests under a surface temper- 
ature varying sinusoidally. 

Excellent agreement is obtained, Figures 1 and 2, between analog results 
and analytic solutions based on the Modified Berggren formula, Equation (2-3), 
for the step change in surface temperature . In Table I, an actual correction co- 
efficient X, equal to the depth of frost penetration X given in the analog solution 
divided by the depth of frost penetration Xg according to the Stefan Equation (2-1), 
is compared with X given by the Modified Berggren formula from Figure 14 ; . 
The agreement is perfect. 

Analog results for tests 3 and 4 under a surface temperature varying 
sinusoidally are compared with results of studies on the IBM-701 computer, 
Table!, since the Modified Berggren formula is not strictly applicable to this 
case. The agreement is not as good as before. However, the errors due to 
"lumping" of the soil profile and especially those caused by the physical per- 
formance of the latent heat wells become larger when the freezing process is 
followed by thawing. Gradients below the frozen zone are very small when thaw- 
ing starts (for example see Figure 3 at a time equal to 250 days) which makes 
the lower boundary of the frozen soil very difficult to determine. 


*The short horizontal lines in Figures 1 through 13 are related to the analog so- 
lution and may be explained as follows. The left end of the line represents the 
time when the soil "lump" begins to freeze, in other words, the time when the 
latent heat well starts to drain. The right end gives the time when the "lump" 
is entirely frozen. In the case of thawing, the reverse is true. 
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3-03. SOLUTIONS FOR STRATIFIED SOIL 

A program of nine solutions on the hydraulic analog computer was car- 
ried out on six different soil profiles, all having a 3-inch bituminous concrete 
surface. Two of these tests were run under a surface temperature varying 
sinusoidally over an annual period. Samples of the scaling and programming 
procedures for analog solutions No. 9 and No. 12 are given in Appendix B. 

Results of the analog solutions are plotted in Figures 5 through 13 t A 
comparison of these results with formulas, summarized in Table II, will be 
made following a discussion of the computational procedures for predicting frost 
depth according to the Stefan and Modified Berggren formulas . 


IV. FORMULAS FOR DEPTH OF FROST PENETRATION 
IN MULTILAYER SOIL PROFILES 


4-01. GENERAL 


No exact analytic solution for the depth of frost penetration in multilayer 
soil profiles exists. However, approximate solutions based on formulas for 
homogeneous soil have been presented. The use of the Stefan equation and Mod- 
ified Berggren formula for this purpose will be presented here. 

4-02. STEP CHANGE IN SURFACE TEMPERATURE 


A step change in surface temperature is defined as that case where the 
temperature of the soil is initially uniform and the surface temperature is sud- 
denly changed from this initial temperature above freezing to some temperature 
below the freezing point. 

a. Stefan Equation. The depth of frost penetration from the Stefan equa- 
tion, which ignores the volumetric heat of the soil, may be written: 


, /48k.v t / 48k-F 
X S = y-lr- = Vtt 

for uniform soil under a step change in surface temperature. For stratified soil 
two approaches yielding identical results may be used: 

(1) "St. Paul" Approach. The so-called "St. Paul" method* makes 
use of the following relationship for the partial freezing index required to freeze 
any layer "n": 



♦This method is presented in Chapter 6, Part XV of the Engineering Manual for 
Military Construction, Corps of Engineers, U. S. Army, dated October, 1954. 
The derivation of Equation (4-2) is also given on Page 44 of Reference 2. 
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where: 

U = L d ; the latent heat removed when freezing the n^ 1 
n n n i a y er of thickness d n « 

d 

R = ; thermal resistance. 

T-l lr 7 


If the freezing index for the locality in question is known, the frost depth can be 
determined by finding those layers whose sum of partial freezing indices equals 
the total for the period. The total freezing index for any period follows from 
Equation (4-2): 



which is Equation (2-54) in Reference 2 when U = L d . 

M v n n n 

the total freezing index may also be written: 


It has been shown that 



(4-4) 


(2) An Equivalent Alternate Approach. An entirely equivalent proce- 
dure to that given above is to determine an effective "L/k" for the frozen zone 
and compute the depth of frost penetration from Equation (4-1). The "St. Paul" 
method is preferred to this procedure because it is more direct. However, this 
alternate approach is fundamental to understanding the nomograph procedure 
(2, pg. 49) where the correction factor is introduced. 

Equation (4-1) for uniform soil may be re-written: 
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We may substitute F from Equation (4-4) for the multilayer case to obtain: 


or: 


Gr) 
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) 
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jn ( m m ) 

k A 2 / 

m N 


(4-5) 


which is the same equation as that given in Step 5, page 49, of Reference 2 . 
Thus, one may estimate the frost depth X, determine the (L/k) g ^ from 
Equation (4-5), then compute X from Equation (4-1). If the computed value of X 
differs appreciably from that originally estimated, the computation may be re- 
peated. 


The procedure outlined above is identical to that suggested for use with 
the frost penetration nomograph (2) except that here, the correction coefficient 
X, has been omitted which is consistent with an approach which ignores volu- 
metric heat. 

b. Latent Heat Plus Voliftnetric Heat of Frozen Soil. As a second some- 
what less approximate approach, the heat removed from each frozen layer as it 
is cooled from the mean annual temperature to 32 F. may be added to the latent 
heat, L R d‘ . More specifically, we assume that the total heat removed from the 
nth layer is now: 


U = L d + C v d 
n n n non 


(4-6) 


Thus, from Equation (4-6) may be used in Equations (4-2) through (4-5) to 
give predicted depths of frost penetration somewhat closer to the actual. How- 
ever, the heat given off while cooling the unfrozen soil and that released by the 
frozen soil as it cools from 32 F. is still neglected. Therefore, the predicted 
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depths are still greater than the actual . 

c. Adaptation of the Modified Berggren Equation. The depth of frost 
penetration from the modified Berggren equation may be written 

X MB= X (4 ' 7) 

for uniform soil. 

In the development of the Modified Berggren formula the surface tempera 

ture was assumed to change suddenly from v degrees above freezing to v de- 

o s 

grees below freezing where it remained constant. This case is referred to as a 
step-change in surface temperature. The relationship among the dimensionless 
parameters ct , p, and X, Figure 14, is strictly applicable only to this case and: 


o C 

a - — and u. = -j- v 
v . L s 

s 

If the soil is stratified, however, and the surface temperature steps 

suddenly from v degrees above freezing to v degrees below freezing, the fol- 
o s 

lowing procedure for determining a frost depth vs time curve is suggested: 

Step 1. Compute weighted values of volumetric heat and latent heat within the 
depth of frost penetration X for which the time is desired: 


C wt = 


L _ = 
wt 


C.d. + C 0 d 0 + + C d 

11 2 2 mm 


L l d l + L 2 d 2 +•••••+ L ra d ni 
X 


(4-8) 


Step 2. Determine a and jx from: 

v 


a = 


o 


wt 

* =r- v s 
Wt 


Step 3. Enter Figure 14, to obtain \. 

Step 4. The time t required to freeze to depth X is now computed from: 

F 


t = 


v \ 
s 
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where F is given by Equation (4-3) or (4-4). Note that this time is equal to 

2 

the time computed from the Stefan equation divided by X . 

Sample computations using this procedure are given in Appendix B. 

4-03. SURFACE TEMPERATURE VARYING SINUSOIDALLY 

In reality a pavement surface experiences daily as well as seasonal 
temperature fluctuations through local variations in air temperature, wind 
velocity, precipitation, solar radiation, etc. Indeed, a sinusoidal variation in 
surface temperature over an annual period more nearly represents the true case. 
Thus, we examine the effect of such a variation on the computational procedures 
for the depth of frost penetration. 

a. Stefan Equation. The depth of frost penetration from the Stefan equa- 
tion is actually independent of the surface temperature variation. Therefore, 
all procedures given in earlier sections are entirely applicable for determining 
the frost depth for any partial freezing index for a freezing season. 

b. Adaptation of Modified Berggren Equation. In computing the maximum 
depth of frost penetration below a typical highway or airfield pavement the follow- 
ing semi -empirical adaptation of the Modified Berggren formula is suggested. 

It is assumed that the pavement freezing index F, duration of freezing period t, 
mean annual air temperature, and the rmal properties of the soil are given. The 
layers making up the profile will be numbered consecutively downward starting 
with the Portland- cement concrete or bituminous concrete layer as No. 1. 

Step 1. Compute an effective ^ from Equation (4-5). 

Step 2. Compute weighted values of C and L within the estimated depth of frost 
penetration from Equation (4-8). 

Step 3. Compute the effective values of cn and m- from: 


a = 


v t 
o 


and 


F = 


C F 
wt 

L 

wt 


Step 4. Determine the correction coefficient X from Figure 14 
Step 5 . Compute the depth of frost penetration from: 
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X = X 

If the computed depth differs appreciably from the assumed depth, 
steps 1 through 5 may be repeated. Seldom are more than two cycles 
necessary. 

A sample computation using this procedure is given in Appendix B. 






V. RESULTS AND CONCLUSIONS 


5-01. DISCUSSION OF RESULTS 

A summary of results of this investigation for the multilayer soil profiles 
is given in Table II. It is believed that the most satisfactory method of compar- 
ing results of these tests is to determine the actual values of the correction co- 
efficient X. defined by: 



where: 

X : "actual" depth of frost penetration at any time (for any freezing 
index) from the hydraulic analog computer 

Xg: depth of frost penetration at the same time (for the same freezing 
index) according to the Stefan formula . 

and compare this value with the ratio: 

X - Xmb 
MBD" Xg 

where: 

Xmb : depth of frost penetration (for the same freezing index) according to 
the Modified Berggren formula.* 

Thus, when equals X- , the Modified Berggren formula yields an "exact” 

result for the depth of frost penetration. 

Several arbitrary values of X for each test have been selected in order 
to make this comparison, the results of which are given in the last column of 
Table II. It can be seen that X^^ is generally smaller than X (negative per- 
centage) which indicates that the Modified Berggren formula is predicting too c 
shallow. More specifically, if the maximum recorded depth X for each of the 
nine tests is selected for comparison, the mean result is - 2.5% with a variation 
of from - 5% to + 1% . 


*The term Xj^g is reserved for the correction coefficient from the Modified 
Berggren formula as given in Figure 14. When the soil profile is uniform, 

X MB = X MBD' 
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5-02. CONCLUSIONS 

Although the number of cases examined in this investigation is limited, 
the writers believe that these tests demonstrate the general validity of the pro- 
posed computational procedure for the depth of frost penetration based on the 
Modified Berggren formula. The adaptation for the multilayer soil profile 

"“V 

proposed herein appears to yield results which are slightly too shallow. Never- 
theless, this error is believed to be insignificant compared to those errors in 
the predicted depth of frost penetration arising from uncertainties in predicting 
the thermal properties of soil, in particular the thermal conductivity, and in 
determining the pavement surface temperature from air temperature and micro- 
meteorological data. 

Therefore, it is believed that future research related to the prediction of 
the depth of frost penetration should be directed toward these problems. 
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Table I 

SUMMARY OF CHECK SOLUTIONS 


A. Step Change in Surface Temperature 


Test 

No. 

Test 

Results 

V 

0 

(deg F. ) 

V 

s 

(deg F. ) 

a 

F 

X _JL 

x s 

X MB 

1 

Figure 1 

12.5 

12.5 

1.0 

0.2 

0.78 + 

0.78 

2 

Figure 2 

1.0 

12.5 

0.08 

0.2 

0.95 

0 95 


B. Sinusoidal Variation of Surface Temperature 


Test 

No. 

Test 

Results 

v q 

(deg F. ) 

V 

s 

(deg F. ) 

a 

M- 

n 

^IBM 

3 

Figure 3 

12.5 

12.5 

1.057 

0.189 

0. 757 

0. 73 ± 

4 

F igure 4 

1 .0 

12.5 

0.078 

0.209 

0.95 

0.96 + 


Soil Profile: 

Homogeneous: y , = 110 lb per cu ft 
u = 9 . 6 % 

Thermal Properties: 

Thermal Conductivity, k = 1.2 Btu per (hr) (ft) (deg F. ) 
Volumetric Heat, C = 24 Btu per (cu ft) (deg F. ) 
Latent Heat, L = 1500 Btu per cu ft 




































Table II 


SUMMARY OF RESULTS OF TESTS ON STRATIFIED SOIL (1) 


A. Step Change in Surface Temperature 


Test 

Soil Profile^ 

Assumed Soil Properties 




' 


x (3) 

, X MB 


X. - \ 

MBD 

No. 


w 


■1 

C 

k 



a 


X 

„ A 


SI 

0 

s 

M- 

MB 

*7MBD Xg 

K ~ X 
X S 

k mbd 

5 

6" Base 

5% 

110 


23 


9.8 

9.0 

1. 09 

0. 221 

2.55 

0.755 

0.84 

0.77 

+ 8% 



7 

120 

120 

1200 

1700 

25 

27 

■|§|t 

DS 




0. 172 

5.07 

0.785 

0.77 

0.78 

- 1% 

6 

"Dry" Subgrade 

10 

1.0 

9.0 

0. 11 

0. 221 
0. 172 

2.55 

5.07 

0.94 

0.95 

0.93 
0. 94 

0.95 

0.99 

- 2% 
- 5% 

mm 

6" Base 

il 

110 

850 

ESI 

1.00 

1.0 

9.0 

0. 11 

0. 166 

3.27 . 

0. 955 

0.96 

0.98 

- 2% 

: mm 

"Dry" Subgrade 

tm 

120 

1700 


1.70 

0. 160 

4.33 

0.956 

0.93 

0.96 

- 3% 


21. 5" Base 

D 

MEM 

; 

25 

wm 

5.0 

10.0 

0.5 

0.242 

2.05 

0.85 

0.87 

0.86 

+ “1& " 

8 

"Dry" Subgrade 

m 



27 

■ 

0. 187 
0. 181 

4.57 

5.63 

0.87 
0. 875 

0. 87 
0. 86 

0. 86 
0.85 

+ 1% 
+ 1% 

m 

6" Base 

21.5" Subbase 
"Wet" Subcrade 

5 

7 

20 

■kg 
w' 3 
■1*9 

EH 

E9 

iipg 

Iwn 


1070 

075~ 

. 0. 221_ 
0. 128 

_2. 55. . 
5.55 

0..855 

0.90 

7 _ —0. 92 . _. 

0.90 

_0, 895 _ 
0.91 

_ .. + 3%, 

1% 

10 

6" Base 

a 

110 

850 

EM 

1. 00 

5.0 

10.0 

0.5 

0. 124 

2.25 

0.90 

0. 89 

0. 98 

- 10% 

"Wet" Suberade 


102 

2900 

mm 

1.70 

0. 109 

3.75 

0.91 

0.90 

0.945 

- 5% 

li 

21. 5" Base 

7 

120 


25 

in 

5.0 



0.242 

2.05 

0.85 

0.92 

0.93 

- 1% 


10.0 

0. 5 

0. 142 

3. 55 

0. 892 

0. 90 

0. 96 

- 6% 


"Wet" Subgrade 

20 

102 

1 

27 





0. 119 

5.60 

0.903 

0.90 

0.94 

- 4% 


B. Sinusoidal Variation in Surface Temperature 


Test 

No. 

Soil Profile ^ 

Assumed 

Soil Properties 

v o 

v s 

a 

. 

M- 

X 

X MB 

\ X MB 

MBD" X g 

x X 

X. - V 

MBD 

w 

v d 

Hi 

C 

k 

x 

s 

V MBD 


6" Base 

5% 

wm 


23 

1.00 










12 

21.5" Subbase 

7 

mm 


25 


5.0 

10.0 

0.5 

0. 126 

5. 80 

0.90 

0.90 

0.94 

- 4% 


"Wet" Subgrade 

20 

E3 

2900 

27 











13 

21. 5" Base 

HV 

wm 

1200 

25 

MM 

5.0 

10. 0 

0. 5 

0. 119 

5.60 

0. 90 

0.90 

0. 91 

- 1% 


"Wet" Subgrade 

. 

wm 

2900 

27 




l_ 








^ Test Results are plotted in Figures bearing the same number as the Test No. 

( 2 ) 

All profiles have 3 inches of bituminous concrete; assume C = 28, k = 0. 80. 
(3) 

v is obtained from Figure 14. 

MB & 


•^a 
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APPENDIX A 

SINUSOIDAL VARIATION IN SURFACE TEMPERATURE 


When the surface freezing index F and the thermal ratio a are known, 
the equivalent sine curve representing the surface temperature variation may be 
computed, although the direct computation involves the solution of a trigono- 
metric equation by trial and error. However, the main part of the computations 
may be carried out once and for all so that solutions to the final equation may be 
obtained directly from a simple chart. 

The sinusoidal surface temperature function (Figure A-l) may be written: 

(A-l) 


v = a sin cot 


where: 


v = temperature variation from the mean annual temperature, deg F . 

a = amplitude of sine curve, deg F. 

“ = m radlans per da y . .. , 

t = time in days 

The deg F. by which the mean temperature is below the freezing point 
during the freezing period is designated by v while the deg F. by which the 
mean annual temperature exceeds the freezing point is v . The thermal ratio a 
is defined as: 


a = 


(A- 2) 


The freezing index F is defined by the following equations: 

-f • t, 

F 


= I (v - v o ) dt 


where (t^ - t Q ) is the duration of the freezing period. Substituting Equation (A-l) 


and noting that v q = a sin cot 0 ., we have: 


A-l 


A- 2 


F = j (a sin cot - a sin u>t_) dt 


which, on integration, becomes: 


F = - (cos co t - cos a) tj).- a sin cot^t^ ' t ) 
Furthermore, the freezing index may be written: 


(A-3) 


v a sinco t 

F = v (t, - t ) = (t, - t ) = — (t, - t ) 

s' 1 o ' 1 o <x 'l o 


(A-4) 


and the following expressions are apparent from Figure A- 1: 


= — - t 

(A-5) 

(Ji 0 

, = - COS 0)t 
1 0 

(A-6) 


Finally, by combining equations (A-3) through (A- 6) and solving for the thermal 


ratio or in terms of cot : 

o 


a = 


( tt - 2co t ) tan cot 
o o 

2 - (tt - 2cot Q ) tan cot Q 


(A- 7) 


For a given value of a the value of t must be found by trial and error. 
However, Figure A- 2 has been prepared such that the value of t may be read' 


directly for a given a . 

Once t is knov 
o 

following equation (A-l): 


Once t is known, the temperature amplitude may be determined from the 


a = 


sincot 


or 


av 


a = 


smeo 


S- = (. “ \ ( l >1 

t \^sincot o y ^ 182.5 - 2 t J 


(A- 8) 


By choosing the time scale for a computer solution, the speed of the 
kymograph, and the temperature scale, .the actual dimensions of the sine curve 
may be determined for programming the solution. 



SINE CURVE FOR ASSUMED ANNUAL VARIATION IN 
SURFACE TEMPERATURE. 



FIGURE A-l 



APPENDIX B 
SAMPLE COMPUTATIONS 


STEP CHANGE IN SURFACE TEMPERATURE 


TEST NO. 9 

Soil Profile and Assumed Thermal Properties: 




w 


L 

C 

k 

3" 

Bituminous Concrete 

-- 

-- 

-- 

28 

0.8 

6” 

Base 

5 

110 

850 

23 

1.0 

21.5" 

Subbase . 

7 

120 

1200 

25 

1.3 


’Wet* Subgrade 

20 

102 

2900 

27 

1.7 


Temperature Conditions: 

v = 5 F 
o 

v = 10 F 
s 

Step change in surface temperature from 37 F. to 22 F . 
Hydraulic Analog Solution: 

The scaling for this problem is given on the accompanying 
'Scaling Sheet' for ’Hydraulic Analog Computer . ’ 

Temperature Scale V g = 0.2 inches per degree F. 

Thermal Energy Scale U = 0. 00195 cu in. per Btu. 

s 

Time Scale T g =0.5 min per day 

Observations made during the solution are shown on the ac- 
companying 'Observation Sheets/ 

The water temperature during the test remained about 0 6 C 
below the desired 28 C, This temperature deviation causes an error of 
about 0.5 percent which is too small to justify a correction. 

Results of the analog solution are plotted in Figure 9. 


B-l 


B-2 


Stefan Equation: 


d 

X = Sd 

k 

L 

R _ d 

R ~k 

S R + f 

n 

2F 

■ 

0.25 

0.25 


0 

0.31 

0.16 

0 

0 

0 

0.50 

0.75 


850 

0.50 

0.56 

10 

10 

1 

1.80 

2.55 

1.3 

1200 

1.38 

1.50 

135 

145 

14.5 

1.50 

4.05 

1.7 

2900 

0.88 

2.63 

447 

592 

59.2 

1. 50 

5.55 

1.7 

2900 

0.88 

3.51 

636 

1228 

122.8 

0.25 

5.80 

1.7 

2900 

0.15 

4.03 

127 

1355 

135.5 

0.42 

6.20 

1.7 

2900 

0.25 

4.23 

215 

1570 

157.0 


* p = n n (zr+I?^ 
n 24 \ + 2) 


See Figure 9 for plotted results. 


Adaptation of Modified Berggren Equation: 


= — = 0.5 
v 

s 


C C ^ 

„ Wt m Wt 

l 1= — v s = 1 ° — 


Wt 


Wt 


d 

X = 2d 

C 

C + 
wt 

L 

L ^ 
wt 


X MB 

2F 



0.25 

0.25 

28 

28 

0 

0 


-- 

0 

0 

0.50 

0.75 

23 

24.7 

850 

567 

0.392 

0.795 

10 

1.6 

1.80 

2.55 

25 

24.9 

1200 

1015 

0.221 

0.855 

145 

20 

1.50 

4.05 

27 

25.7 

2900 

1712 

0.150 

0.890 

592 

75 

1.50 

5.55 

27 

26.0 

2900 

2030 

0.128 

0.897 

1228 

153 

0.25 

5.80 

27 

26.1 

2900 

2070 

0.126 

0.90 

1355 

167 


See Figure 9 for plotted results. 
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SINUSOIDAL VARIATION IN SURFACE TEMPERATURE 


TEST NO, 12 

Soil Profile and Assumed Thermal Properties: 

(Same as Test No, 9) 

Temperature Conditions: 

v =5F 
o 

v = 10 F 
s 

Sinusoidal variation in surface temperature 


Surface Temperature Curve: 
From Figure A- 2: 


v 

a - — - 0.5 therefore t =14,35 days 
v o J 


F = v g (182,5 - 2t Q ) = 1538 degree days 


a = 


o 


sin cot . 2-rr , . oe 
o sm^j 14.35 


= 20.45 F 


Therefore, the following equation is obtained for the surface temperature 
variation: 


v = 20,45 sin t 
365 


Hydraulic Analog Solution: 

The scaling for this problem is the same as that for Test 
No. 9 and is therefore given on the accompanying ’Scaling Sheet' for 
Test 9 . 

Observations made during the analog solution are shown on 
the accompanying ’Observation Sheets'. Data were taken during one 
surface temperature cycle . The initial temperature distribution in the 



ground was assumed on the basis of results of Test No. 3 which repre- 
sents uniform soil. This apparently crude approximation was later found 
to yield errors which are very small. 

It was difficult to get a good determination of the maximum 
depth of frost penetration in this test. Shortly after the thawing started 
at ground surface, gradients were very small at the freezing interface . 
Latent heat well No. 11 was 90 per cent drained and well No. 12 had just 
started to drain. As thawing continued downward from the surface, well 
No. 12 continued to drain into well No. 11 which complicated the usual 
procedure of determining the maximum depth of frost penetration. How- 
ever, it is felt that the most reliable determination of the maximum depth 
of freezing can be obtained from observations at the start of the thawing 
period when 70 per cent of well No. 11 was frozen. Thus: 

X = (5.61 - 0.51) + (0.7) (1.02) = 5.81 ft . 

The actual surface temperature curve as observed during the 
test deviated slightly from that intended yielding the following results: 

F = 1568 degree days 

t = 157 . 5 days 

Therefore, actually: 

v = 9. 95 F 
s 

S' = 0.503 

fj. = 0.126 (see computations for Test No. 9) 

Results of the analog solution are plotted in Figure 12. 

Stefan Equation: 

From sample computations for Test No. 9, the depth of frost 
penetration corresponding to F = 1568 degree days is: 

X g = 6.20 ft 

Therefore: 



B-5 



5.81 

6.20 


0.94 


A daptation of Modified Berggren Equation: 


We start with Step 1 of the procedure suggested in Section 


4-03-b: 


Step 1. 

Assume that the depth of frost penetration is 5. 7 feet. Thus: 



2 

(5.7) 2 


■ n 425 2150 9160 

(0 + (850) (0.5) + (1200) (1 . 79) + (2900) (3. 16) ) 


+ (212 +2150 + 9160) 

+ T1§ (1075 + 9160) +j^ (4580) 



= 1970 


Step 2. 


C = (28) (0.25) + (23) (0.5) + (25) (1.79) + (27) (3.16) = 26 0 
wt 5.7 


L = 0 + (850) (0.5) + (1200) (1. 79) -f (2900) + (3. 16) = 2Q6Q 

wt 5.7 


Step 3. 


OL 


V (5) (157.5) 
F 1568 


0.503 


M- = 


C 

wt 

L J 

wt 


(26.0) (1568) 
(2060) (157. 5) 


= 0.126 


Step 4. 


From Figure 14, = 0.90 
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HYDRAULIC ANALOGUE COMPUTER 
SCALING SHEET 

PROBLEM 1 4" layer Case, (Vef Subgrade 


TEST NO.: *7 

DATE : Jo/ y 20. rfS5 


COMPUTED BY: (2/Sf 



f2.£3is-h/' Type A H-* ■Xee/J/or' Type 3 































































HYDRAULIC ANALOGUE COMPUTER observed by: test no ^ 

4 - layer Case - (Vet SubgraVe OBSERVATION SHEET date = Jvty2o,i4ss sheet L- of ^ 

PROBLEM •' Step Change /b> ^Surface Te/r?per a tore 


DATE 

AND 

HOUR 

ELAPSED TIME 
JLmIN. r^-OAYS 

STANDPIPES : TEMPERATURE SCALE a.2 INCH = 1° F 

REMARKS 
•*- Depth tn Ft. 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

m 

12 

m 

IB 

IB 

IB 

17 

18 

19- 

20 

ANA- 

LOGUE 

(MIN) 

PROTO 

TYPE 

(DAYS) 

HH 

0- /3 

0.36 

0.G3 

04 7 

142 

iec> 

2-31 

24& 

532 


SGI 

Ml 

14 Z 

f.So 

//.46 

/£?4 

IB 

BIB 

5022 

6/42 

13- Za 

O 

O 

2 Z 

31 

31 

m 

ESI 

m 

31 

31 

31 

eh 

37 

m 

m 

31 

31 

HI 

*7 

3? 

EH 

37 

^2 







s" 

3* 

hoo 

2-35 

4-4o 

IB 

!o35 

H ss 

32 os 

BB 

m 

— 

su 



IB 


b -fret 








25" 

So" 

2-3o 

4-4o 

TM 

ll-IO 

EB 


(,8:3c 

to-Z- 

jzz-do 


HI 


m 

* 

PjB 

■Jy -fr 

oze»l 



n 

to 

zz 

225 

2f. 5 

26." 



3Z 

325 

33 



3*1 

35.2 

355 

36-Z 

36.7 

37 .- 

HI 

HI 

57 

?7 



to 

Z» 

ZZ 

Qi 

24.1 

25.2 

26.7 

4?.? 

m 

312 

32- 

325 

13 

134 

33.1 

345 

35. Z 

362 

56.7 

365 

EH 

37 

37 



i5 

3o 

ZZ 

tz.t 



26-2 

27.6 

22.0 


3/7 

32.2 

32-7 

33' 

J3.5 

346 

346 

357 

366 

36? 

37 

37 

T7 

714 >°c CTcSt tc*tp) 


Za 

4-0 


224 

23.5 

Z44 

25.7 

*7 

28.3 


3/ / 

32 

$Z.4 

377 

33.1 

33.7 

34/ 

335 

36* 

36.7 

37 

37 

37 



25 

So 

ZZ 

224 

234 

24-3 

25.5 

26.6 

274 

m 

342 

32 

3*3 

32-6 

32-7 

33.4 

33.1 

35 

56 

365 

BE 

37 

37 

17 . 3 °4 


3o 

6 >o 

ZZ 

22 S 

234 

24. Z 

25.4 

26.6 

a 7-7 

3g/j 

3a 1 

3/.7 

322 

325 

BH 

332 

337 

344 

357 

364 

Bl 

57 

3? 



35 

lo 

zz 

2ZS 

23.4 

24 . 1 

2S.2 

23,3 

274 

27.5 

201 

3/-4 

32 Z 

3z5 

3Z.6 

33.1 

33.6 

345 

35.6 

36* 

36.1 

? 7 

»? 

Z7.4C 


4o 

So 

ZZ 

Z24 

23.4 

242 

252 

74.3 

214 

> 8 . 5 

100, 

3/. 2 

32.2 

325 

32-1 

33.1 

355 

BH 

35 3 

36.' 

36. f 


37 



4-5 

lo 

zz 

zze 

23.5 

24 3 

252 

2(r.Z 

2 7-2 

27.3 

so 

HI 

322 

32-3 

3*7 

336 

335 

?4.2 

3S-I 

35? 

36.? 

361 

37 



So 

too 

zz 

zz-5 

*3.2 

24 


254 

SB 

RE 

227 

30.4 

3/7 

3 2- 

323 

327 

35. 1 

54/ 

35.1 

35.7 

36.7 

56 f 

37 

274 °C 


55 

no 

zz 

tZ-3 

23 

234 

24-7 

25.6 

26.7 

27.6 

27.6 

3a 

3/5 

3v 

32.2 

32-7 

33 

56 

341 

356 

367 

364 

57 



(?o 

tzo 

zz 

Z2-1 

23.5 

242 

243 


Z 67 


28 . 5 

27 . f 

31- * 

HU 

37.4 


33 

33? 

34? 

354 

36.7 

161 

37 



ss 

Ho 

_£^ 

2Z3 

23. Z 

23. S 

246 

25 5 

264 

274 

27 3 

BS 

3/. 2 

32.2 

31 4 

32 7 

33 

335 

347 

35 3 

367 

56.7 

37 



lo 

<4o 

zz 

23 

23-5 

Z4-i 

247 

256 

26.5 

27.3 

27.2 

EH 

30.7 

32 

524 

3z-6 

ESI 

337 

34.6 

35* 

365 

367 

V 



IS 

/SO 

zz 

224 

23.1 

23.5 

2f 5 

257 

262 

H 

27.7 

27 z 

3 o. S 

31-7 

32-1 

325 

3?5 

EH 

345 

35./ 

^64 

361 

57 



0 O 

/ 6 >o 

zz 

223 

23 

23.5 

245 

253 

*62 

27 

27? 

*7 

30 3 

31-6 

3Z 

>2-3 

3*7 

33.6 

544 

35 . 1 

364 

367 

57 

21S e C 


es 

4° 

zz 

22 . S 

23 .1 

23 5 

244 

25.2 

26 

26.7 

277 

t%4 

30.1 

31-3 

5Z 

323 

32-6 

35.5 

343 

341 

363 

361 

37 



0 o 

> to 

zz 

224 

23 

23.5 

2f.3 

zs 

252 

26. 7 

27.5 

2 fl 

247 

3l-> 

3 Z 

323 

32.6 

33.4 

54* 

347 

362 

36i 

>L_ 



is 

ifo 

zz 

223 

22- ‘ 4 

23.5 

243 

25 

* 5 ? 

26-7 

27.5 

27.6 

24? 

3/1 

32 z 

32.4 

JZ.6 

335 

341 

347 

36> 

36? 

57 



/oo 

Zoo 

zz 

22 S 

2.3 - 1 

23.5 

2f.3 

25/ 

25.? 

267 

274 

28 . S 

20. (0 

30 i 

3Z 

3Z3 

32.6 

333 

34 » 

34.7 

36-2 

36.? 

37 



loS 

Zlo 

zz 

Zz3 

23 

23.5 

244 

25. 1 

25. S' 

26-7 

274 

275 


5>.l 

EH 

32 2 

325 

53 .Z 

34' 

546 

36.i 

36. S 

37 
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